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A momentum-space approach of the density-matrix renormalization-group (DMRG) method is de- 
veloped. Ground state energies of the Hubbard model are evaluated using this method and compared 
with exact diagonalization as well as quantum Monte-Carlo results. It is shown that the momentum- 
space DMRG is a very useful numerical tool for studying the Hubbard model and other fundamental 
models of interacting electrons in two dimensions. For the Hubbard model in two dimensions, the 
momentum-space DMRG method reproduces accurately the exact diagonalization results of ground 
state energies on a 4 x 4 lattice and yields new upper bounds of ground state energies on an 8x8 
lattice. 



PACS number: 02.70.Rw, 71.27.+a 

Thaf eal-space DMRG method, which was proposed by 
Whitetl in 1992, is a powerful method for studying ground 
state properties of interacting electrons or spins. It has 
been successfully applied to various one-dimensional (IIX) 
quantum models, such aSr4he spin Heisenberg modelM, 
the Kondo lattice modelcl, the Hubbard modelcl, and 
other models!! n However, in 2D its application is still 
not satisfactory^. In this paper we generalize the DMRG 
method to momentum space. We hope this may enlarge 
the range of the application of the DMRG method and 
provide an accurate numerical method for studying inter- 
acting electrons in 2D. We shall take the Hubbard model 
as an example to show how the DMRG method works in 
momentum space. The properties of the Hubbard model 
have been extensively studied in both one and two dimen- 
sions, and so this provides a good background for com- 
paring the momentum-space DMRG method with other 
methods. 

The DMRG method is a diagonalization technique 
which attempts to use a small number of states, say m 
states, to expand the ground state (or some low energy 
excitation states) accurately. In conventional numerical 
renormalization-group (RG) methods, one keeps_the m 
lowest energy eigenstates of a block Hamiltonianllia. In 
the DMRG method, however, one keeps the m most prob- 
able basis states in describing the ground state of a larger 
block, called a superblock. A superblock contains two 
blocks, a system block and an environment block. In a 
DMRG iteration, one diagonalizes the Hamiltonian of the 
superblock, finds out the reduced density matrix for the 
system block from the ground state of the superblock, 
truncates the Hilbert subspace by keeping the m largest 
eigenstates of the density matrix, adds one or more sites 
to the system block to form a new superblock, and then 
repeats the above procedure until the desired result is 
obtained. 

There are two approaches in constructing a superblock, 
an infinite-lattice approach and a finite-size approach. 
In the infinite-lattice approach, the environment block 
is generally chosen as the space reflection of the system 
block. In the finite-size approach, on the other hand, the 



size of the superblock is fixed and the environment block 
is chosen as the remaining part of the lattice for a given 
system block. The infinite-lattice approach allows the 
size of the superblock to be flexible and can be used to 
study directly the thermodynamic limit. The finite-size 
approach is, however, more accurate in calculating quan- 
tities for a finite lattice system. In White's original paper, 
two sites are added to a superblock each time, one to the 
system block and the other to the environment block. 
This way of constructing a superblock preserves the re- 
flection symmetry and avoids some problems caused by 
odd-size lattices if the infinite-lattice approach is used in 
ID. However, it is generally more efficient in computation 
if one adds just one site to the system block and no site 
to the environment block when the finite-size approach 
is used. 

The DMRG method satisfies the variational principle 
because it uses a small physical subspace to approximate 
the full Hilbert space and no unphysical states enter the 
truncated Hilbert subspace. Thus the ground state en- 
ergy obtained using this method is always a upper bound 
for the true value. The error for the ground state energy 
decreases as the number of retained states is increased. 
The truncation error is generally smaller than the true 
error of the result. 

The Hubbard model is defined by the Hamiltonian 
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where < ij > means summation over nearest neighbors. 
In momentum space, it reads 
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with TV the lattice size. Here the periodic boundary con- 
dition is assumed. Each momentum k point has four de- 
grees of freedom, i.e. { |0 >, c\^\0 >, c\^\0 >, c^c^JO > 
}. In the RG iteration, generally a k point with these 
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fours states can be added to a system block to form a 
new system block. However, in our calculation we treat 
the spin degree of freedom as an extra spatial coordinate 
and a momentum-spin point (fee), which has only two 
degrees of freedom {|0 >, c\ a \0 >}, as a basic unit which 
is added to a system block. In this case, the spin rotation 
symmetry is broken, but the total number of degrees of 
freedom of a superblock is reduced and more states can 
be retained at each truncation of Hilbert subspace. This 
can reduce truncation errors and save computer time. 

For the convenience in the discussion below, we label a 
superblock as (A • B), where • represents a momentum- 
spin point (kcr) which has just been added to the system 
block, A represents the collection of all the momentum- 
spin points in the system block excluding the point (fcer), 
and B represents the collection of all the momentum-spin 
points in the environment block. There are many ways in 
ordering k points in a system. In calculation, one should 
test a few possibilities of orders of k points and choose 
the one which gives the lowest ground state energy. 

The Hubbard interaction is local in real space. How- 
ever, it is non-local in momentum space; it contains terms 
which link two or more k points in momentum space. The 
summation in the second term of (Q) contains iV 3 terms. 
They can be grouped into the terms which are defined 
purely in each sub-block and the terms containing inter- 
actions among A, • , and B. To evaluate matrix elements 
of the Hamiltonian, we find it is very useful to define the 
following composite operators of electrons: in block A, 
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where a = —a; in block B, the corresponding operators 
can be obtained by changing (a, A) into (b, B). They are 
the basic operators whose matrix elements are kept and 
updated in our DMRG iterations. The total number of 
these composite operators is 6N at each sub-block. Us- 
ing these composite operators, the Hamiltonian (|^) can 
be rewritten as H = H, + Ha + Hb + Ha, + Hb, + 
Hab + Ha»b, where H % = ekUka-, Ha is the same as the 
Hamiltonian (^) with Cka replaced by ao(ka), 

H A . = — jn fe(T ai(0,o-) + 4 CT a 3 (fco-) + a\{ka)c ka ^ , 
Hab = Jf E{E [^(pa'M-pa') + bl(j>a')a 3 (pa') 

p a' 

+(a < — > b)] + a\{p)bi{p) - b 2 (p)4(p)\ + h.c, 



Ha.b = jjc\ a 2^|ai(fc -pa)b (pa) - a 2 (p - k)b {po) 
p 

+abl(pa)a4,(k + p)) + (a < ► b)\ + h.c, 

and Hb and Hb* can be obtained from Ha and Ha% by 
changing a into b. In Ha»b, o- 2 (p) = a 2 (p) if f =T or 
a 2 {p) if o =1. 

The breaking of lattice symmetries in truncating ba- 
sis states is a unavoidable source of errors in the DMRG 
method. It is also the main difficulty in using the DMRG 
method to in real space in 2D. In real space, the transla- 
tion symmetry, i.e. the momentum conservation, is bro- 
ken in the DMRG method. In momentum space, how- 
ever, this symmetry is preserved. This symmetry prop- 
erty is undoubtedly very useful. Combined with two 
other conserved quantities, the number of up spins iVj 
and the number of down spins Ni, it can be used to 
block diagonalize the Hamiltonian. Not only can this 
save the computer time, but also it allows us to keep 
many more eigenstates in the truncation of Hilbert space 
as the number of non-zero matrix elements is now signif- 
icantly reduced. Basis states at each block can be clas- 
sified by three quantum numbers (iVj , JVj, , P) with P the 
total momentum. For the composite operators defined 
in (||) it can be shown that their matrix elements are 
nonzero only when the difference between the momen- 
tum of the initial state and that of the final state is p. 

In the DMRG method, the basis states for both the sys- 
tem block A and the environment block B are incomplete. 
If the Hamiltonian contains terms with interactions be- 
tween A and B, the matrix elements of these terms will 
be less accurately approximated after the truncation of 
the Hilbert space compared with the matrix elements of 
other terms which arc defined purely within each block. 
In ID in real space one can choose A and B so that no in- 
teractions exist between A and B, i.e. Hab = Ha»b = 0, 
if all terms in the Hamiltonian arc local in space. In real 
space in 2D or in momentum space in any dimension, 
however, interaction terms between A and B always ex- 
ist no matter how A and B are constructed. Thus in 
general results obtained by the DMRG method in real 
space in 2D or in momentum space in any dimension will 
not be as accurate as, for example, the ground state en- 
ergy of the ID spin-1 Heisenberg model that White, and 
Huse obtained using the real-space DMRG methodnl. 

In momentum space, different size lattices have dif- 
ferent k points. Thus in using the DMRG method in 
momentum space, the size of the lattice needs to be fixed 
at the beginning. This means that only the finite-size 
approach of the DMRG method can be used in momen- 
tum space. To use the finite-size approach, however, one 
needs first to build up a series of initial system blocks 
and the corresponding environment blocks. As the infi- 
nite lattice approach of the DMRG is not applicable, in 
this case, we shall use the conventional RG methodEm to 
build up the initial system and environment blocks. The 
steps in building up the initial system blocks (similarly 
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for the environment blocks) are as follows: 1. Start from 
a small system block A±, which can be handled with- 
out truncation of basis states. 2. Add a new (fcer) point 
to A\ to form a new system block A2. 3. Diagonalizc 
the Hamiltonian in the Hilbert space spaned by A2. 4. 
Truncate the Hilbert space by retaining m lowest energy 
eigenstates, but restrict the number of states retained at 
each (iVf , N±, P) subspace not more than a small integer 
n. Here n is a variational parameter which should be 
determined so that the final result for the ground state 
energy is minimized. In our calculations, we find that 
n = 1 or 2 generally gives the best result for the ground 
state energy for 8 x 8 systems with 1000 states retained. 
The reason for limiting the number of retained states 
at each (JVf,iVj,,P) subspace is to prevent the retained 
states from being centralized in a few (iVf , 2Vj_, P) sub- 
spaces in A2 . Otherwise, some of the (iVf , JVj., P) sub- 
spaces in A2, which may have substantial contribution 
to the final ground state energy, may be neglected in the 
DMRG iterations later. 5. Replace A2 by A\ and repeat 
the steps 2-5 until all the initial system blocks required 
are established. 

All (iVf ,N±,P) subspaces that a system block can 
have are determined purely by the (ko~) points the block 
contains. In the above initialization step, in fact not 
all (-/Vj, Ni, P) subspaces in a system (or environment) 
block need be considered. If a subspace with the quan- 
tum numbers (Ni^,Ni^,Pi) in a system (environment) 
block can not find a subspace with the quantum numbers 
(A^2,T) -^2,1) P2) in the corresponding environment (sys- 
tem) block such that (iVi >T + N 2 ,t,Ni tl + N 2 ,i,Pi + P 2 ) 
is equal to the quantum numbers of the ground state 
required, then this (Ni^,Ni t i,Pi) subspace in the sys- 
tem (environment) block will have no contribution to the 
ground state and can be ignored. 

After the above initialization step, the finite-size ap- 
proach of the DMRG method will be used to find out the 
eigenvalue and eigenfunction of the ground state for a 
given filling factor and momentum. Now m largest eigen- 
states of the density matrix are retained at each trunca- 
tion of Hilbert space and the number of retained states 
at each (N^N^, P) subspace will no longer be limited. 

We have evaluated ground state energies of the Hub- 
bard model in both ID and 2D using this momentum- 
space DMRG method. In most of our calculations 1000 
states are retained at each truncation of basis states. The 
number of states retained at each [N-\ , JVj_ , P) subspace 
on average is small. Our results are much better than 
those obtained by the conventional RG methodB. For a 
10-electron system with U=4 (the energy is measured in 
unit of t, i.e. t=l) in 2D, for example, the momentum- 
space DMRG result for the ground state energy with 1000 
states kept is —19.57, which is much lower than the coru 
ventional RG result with even 3000 states kept, — 18.541B 
(their relative-errors with respect to the exact diagonal- 
ization resultEJJ are 0.05% and 5%, respectively). 

Table I compares the ground state energies obtained 
by the DMRG method in momentum space with those 



obtained in real spaced at half-filling in ID. The ground 
state energies obtained by the DMRG method are lower 
in real space than in momentum space if the same number 
of states is kept, which means that in ID the DMRG 
method works better for the Hubbard model in real space 
than in momentum space. This is not surprising because 
the Hubbard interaction is a local interaction and at half 
filling all electrons are localized in space as a result of a 
Mott insulator transition. 

Table II compares the momentum-space DMRG results 
with the exact diagonalization and the quantum Monte- 
Carlo results on 2D square lattices. The largest lattice we 
have studied so far is 12x12. Compared with the exact 
results on 4x4 lattice, we find that the DMRG results 
are very accurate when U is small. The relative error 
for the DMRG result with U = 2 is 3 x 10" 4 . The mo- 
mentum space DMRG method works better in the weak 
coupling limit because the single particle basis state used 
here is the plane wave state. When U=0, the ground 
state is a filled Fermi sea of non-interacting electrons, 
the momentum-space DMRG method gives the exact re- 
sult for the ground state even when only one state is kept. 
For large U, the DMRG results are not as good as in the 
weak coupling limit. However, they are still comparable 
with those obtained by the projected quantum Monte 
Carlo and the stochastic diagonalization methods. Both 
the DMRG method and the stochastic diagonalization 
method satisfy the variational principle (the projected 
quantum Monte Carlo method does not satisfy the vari- 
ational principle due to the important sampling), and 
the results obtained from these two methods give upper 
bounds of ground state energies. On a 8 x 8 lattice, the 
ground state energies obtained by the DMRG method 
are systematically lower (hence better) than the stochas- 
tic diagonalization results. The DMRG results therefore 
provide new upper bound for the ground state energies 
of the Hubbard model on these lattices. 

In ID, no exact diagonalization results are available 
for the 16-site Hubbard model, but the results obtained 
with the real space DMRG method are very accurateO. 
If we use the best result for the ground state energy of 
the ID 16 site Hubbard model obtained by the real-space 
DMRG method as the true value of the ground state en- 
ergy, we estimate that the relative error for the ground 
state energy obtained by the momentum-space DMRG 
method with m=1000 is 5% when U=4, which is higher 
than the corresponding value in 2D. Thus the momen- 
tum space DMRG method works better for the Hubbard 
model in 2D than in ID. Physically this is because the 
contribution of the kinetic energy term to the ground 
state, which is rigorously treated in the momentum-space 
DMRG method, is larger in 2D than in ID. 

The above discussion for the momentum-space DMRG 
method has been focused mainly on the Hubbard model. 
However, it can be easily generalized to apply to sev- 
eral other physically interesting models, such as the 
Anderson lattice model and the interacting fermion 
model with nearest-neighbor Coulomb potentials. For 
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the Anderson lattice model, the composite operators 
defined in (|§|) can be used without modification. In 
other cases we need to generalize the definitions of 
the composite operators in (|j|) (for example 03 in 
(§) should be defined as a 3 (pa) = Y. qiq2 a> V (<H ~ 

q2)al(qi<J l )a (q 2 (j')ci Q (p + q 1 -q2<T)6(p a)? :A for a gen- 
eral electron-electron interacting model ^2 V{Q ~ 
q')c qtT c q ' cr c i q „ +q , cr ,c q ii +qa /) and introduce some new com- 
posite operators of electrons. For an arbitrary inter- 
acting fermion model with a Coulomb-type potential, 
the number of composite operators required is gener- 
ally of order N 2 , which will limit the application of 
the momentum-space DMRG method to small lattices. 
However, if V(q — q') can be factorized as a sum of 
products of a function of q and a function of q' , i.e. 
V{q~ q') = Ya f 1 (q) 9 1 {<!')■> ^ can be shown that the total 
number of composite operators required can be reduced 
to the order of the system size N. In that case a 12 x 12 or 
even larger lattice system is accessible by the momentum- 
space DMRG method with presently available computer 
facilities. The Hubbard model is obviously a factorizable 
potential, with fi = g\ = 1 and n = 1. The nearest- 
neighbor Coulomb potential is also a factorizable poten- 
tial: V(q — q') ~ cos(g — q') = cos q cos q' + sing sing'. 
For this nearest-neighbor Coulomb potential, the total 
number of composite operators needed is 5N in ID and 
7N in 2D if fermions are spinless. 

In conclusion, we have generalized successfully the 
DMRG method to momentum space and studied ground 
state properties of the Hubbard model. Our results show 
that the momentum-space DMRG method is a power- 
ful numerical method for studying the Hubbard model 
in 2D. There is no reason to believe that this is the only 
model for which the momentum-space DMRG method 
can work. For any finite ranged potential, if it is fac- 
torizable, we believe that the momentum-space DMRG 
should work even better than for the Hubbard model. 

The numerical calculations reported here were per- 
formed on a HP735/99 workstation. In the DMRG it- 
eration, all intermediate data were stored in a hard disc. 
For the calculation on an 8 x 8 (12 x 12) lattice with 1000 
states kept the computer memory space needed is 36 (45) 
Mb in RAM and 200 (500) Mb in hard disc. 

We would like to thank G. A. Gehring and R. J. Bursill 
for fruitful discussions and Xiaoqun Wang for supplying 
the real-space DMRG results in Table I. 
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TABLE I. Comparison of the ground state energy per site 
obtained using the DMRG method in momentum space with 
that obtained in real space at half-filling on a 16-site chain. 





Momentum 


space 


Real 


m 


U=l 


U=4 


U=l 


400 


-1.02925 


-0.51316 


-1.02958 


600 


-1.02944 


-0.53574 


-1.02969 


800 


-1.02952 


-0.53724 


-1.02972 


1000 


-1.02958 


-0.54562 




1200 


-1.02959 


-0.55218 





U=4 



-0.575896 
-0.575900 
-0.575901 



TABLE II. Comparison of the ground state energy ob- 
tained using the DMRG method with m=1000 with the exact 
results on 4x4 lattices (Exact), the cluster diagonalization 
results on 6x6 lattices (CD), the projected quantum Monte 
Carlo (QMC) and the stochastic diagonalization (SQ) results 
on 2D square lattices with N electrons. 



L X X Ly 


U 


N 


ExactEa/CDE 


DMRG 


QMC 


La 


sqN 


4x4 


2 


16 


-18.01757 


-18.012 








4x4 


4 


14 


-15.74459 


-15.673 








4x4 


4 


16 


-13.62185 


-13.571 


-13.6 


-13.59 


4x4 


8 


16 


-8.46887 


-8.263 


-8.4* 


\ 




6x6 


4 


26 


41.49 


41.108 


41.9 


1 


40.77 


8x8 


4 


10 




-34.325 


-34.3 


-34.31 


8x8 


4 


18 




-54.394 


-54.6 


-54.37 


8x8 


4 


26 




-66.098 


-66.8 


-66.05 


12x12 


4 


18 




-64.107 
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